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ABSTRACT 

We present the results from a sensitive multi-wavelength analysis of the properties of Extremely Red Objects 
(EROs). Our analysis employs deep RIzJHK photometry of a 8.5' x 8.5' region to select a sample of 68 EROs with 
(R— K) > 5.3 and brighter than K = 20.5 (5ct). We combine this photometric dataset with an extremely deep 1.4- 
GHz radio map of the field obtained from the VLA. This map reaches a 1-er limiting flux density of 3.5/xJy making 
it the deepest 1.4-GHz map taken and is sensitive enough to detect an active galaxy with L1.4 > 10 23 WHz" 1 at 
z > 1. If powered by a starburst, this radio luminosity is equivalent to a star-formation rate of > 25 M Q yr" 1 for 
stars more massive than 5M Q . We identify radio counterparts to 21 of the EROs in this field with radio fluxes 
above 12.6/iJy and resolve a third of these with our 1.6" FWHM beam. The spectral energy distributions of 
the majority of these galaxies are consistent with those expected for dust-reddened starbursts at z ~ 1 . At these 
redshifts the radio luminosities of these galaxies indicate a median far-infrared luminosity of this population of 
Lfir Z 1O 12 L0, meaning half of the radio-detected sample are ultraluminous infrared galaxies (ULIRGs). We 
conclude that > 16 ± 5% of the ERO population brighter than K = 20.5 are luminous infrared galaxies (LIRGs) at 
z ~ 1 . We also use photometric classification of the colors of the EROs to investigate the mix of dusty active and 
evolved passive systems in the remaining ERO population which is undetected in our radio map. Based on this 
we suggest that at least 30% and possibly up to ~ 60% of all EROs with (R— K) > 5.3 and K < 20.5 are dusty, 
star-forming systems at z Z 1 ■ Our best estimate of the star formation density in this highly-obscured and optically 
faint (R > 26) population is (0. 1-100 M Q ) = 0.1 1 ±0.03 M yr" 1 Mpc" 3 , comparable to estimates of that in Ha 
emitting galaxies at z ~ 1, and greater than the estimates from UV-selected samples at these epochs. This lends 
support to the claims of a strong increase in the contribution from obscured systems to the star formation density 
at high redshifts. Using the observed counts of the radio-detected ERO population we model the apparent break 
in the /^-band number counts of the whole ERO population at K ~ 19-20 and propose that the passive ERO class 
dominates the total population in a relatively narrow magnitude range around K < 20, with dusty, active EROs 
making up the bulk of the population at fainter limits. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: starburst — infrared: galaxies 
galaxies: clusters: individual (Abell 851, CI 0939+47 13) 



1. INTRODUCTION 

The last five years have seen a growing appreciation of the 
diversity of galaxy properties at z Z 1-3. In part this has arisen 
from an impressive improvement in the quality of the observa- 
tions of galaxies at these redshifts - driven by the availability 
of powerful new instruments on 8-m class telescopes (e.g. NIR- 
SPEC on Keck and ISAAC on the VLT). However, an equal role 
has been played by the realization of the necessity of a multi- 
wavelength approach to studies of galaxy evolution at z > 1. 
Thus a range of new surveys spanning wavebands from the X- 
ray, near- and mid-infrared, submillimeter and out to the ra- 
dio, have complimented the traditional view of distant galaxies 
based on UV/optical observations. These new studies have all 
tended to stress the role of dust obscuration in censoring our 
view of the galaxy population at high redshifts, and especially 
in disguising the extent of activity in the most luminous sys- 
tems, both AGN and star-forming. Indeed, our current, incom- 
plete, knowledge of the evolution of dusty galaxies suggests 



that some of the most active galaxies in the distant Universe 
may be so obscured as to be undetectable in even the most sen- 
sitive UV/optical observations (e.g. Haarsma et al. 2000; Smail 
et al. 2002). 

One central theme to appear from this new multi-wavelength 
view of galaxy formation is the ubiquity of optically-faint, but 
bright near-infrared, counterparts to sources identified in many 
wavebands: in the X-ray (Alexander et al. 2001, 2002; Cowie 
et al. 2001; Page et al. 2002; Stevens et al. 2002; Mainieri et al. 
2002), the mid-infrared (Pierre et al. 2001; Smith et al. 2001; 
Franceschini et al. 2002), submillimeter (Smail et al. 1999a; 
Lutz et al. 2001; Ivison et al. 2002) and radio bands (Richards 
et al. 1999; Chapman et al. 2001). This has resulted in a re- 
naissance in interest in such Extremely Red Objects (EROs) - 
a class of galaxies which had previously been viewed as a cu- 
riosity with little relevance to our understanding of galaxy for- 
mation and evolution, comprising as they do a mere ~ 5-10% 
of the population at K < 20. 
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The Properties of ,uJy EROs 



The ERO class is photometrically-defined - the most fre- 
quently used definition is now (R-K) > 5.3. This very red 
optical-near-infrared color is intended to isolate two broad pop- 
ulations of galaxies at z Z 1 : those which are red by virtue of 
the presence of large amounts of dust (resulting from active star 
formation), as well as passive systems whose red colors arise 
from the dominance of old, evolved stars in their stellar popu- 
lations. The very different natures of these two sub-classes has 
prompted efforts to disentangle their relative contributions to 
the ERO population (Pozzetti & Mannucci 2000; Mannucci et 
al. 2002) so that well-defined samples can be used to test galaxy 
formation models (Daddi et al. 2000; Smith et al. 2002a; Firth 
et al. 2002). 

Studies of the ERO population have turned up hints of the 
diverse nature of this population, including an apparent break 
in the number counts of EROs at K ~ 19-20 (McCarthy et al. 
2001 ; Smith et al. 2002a). This feature is roughly 2 mags fainter 
than the well known break in the total TT-band counts of the faint 
field population (Gardner et al. 1993). The cumulative count 
slope of EROs (log 10 JV(< K) oc aK) drops from a = 0.81 ±0.05 
at bright magnitudes to a = 0.32 ± 0.05 for the fainter EROs 
(e.g. Smith et al. 2002a; Chen et al. 2002; Fig. 2). This break 
in the counts could either result from a relatively narrow red- 
shift range probed by the selection criteria of EROs - with the 
change in the integrated counts then being simply related to the 
form of the luminosity function of the ERO population, or it 
may reflect differing evolutionary behavior in sub-populations 
- where one class of EROs declines rapidly beyond K ~ 19. To 
test this suggestion we need to track the relative mix of galaxy 
types within the ERO population to K > 19-20. 

Techniques to separate the passive, evolved and obscured, 
active populations include near-infrared photometric classifica- 
tions (Pozzetti & Mannucci 2000; Smith et al. 2002a; Mannucci 
et al. 2002), submillimeter emission (Dey et al. 1999; Mohan et 
al. 2002), X-ray emission (Alexander et al. 2001, 2002; Brusa 
et al. 2002) and optical/near-infrared spectroscopy (Cimatti et 
al. 2002; Smith et al. 2001, 2002b). Of these approaches, by 
far the largest and most influential project is the "K20" survey 
undertaken by Cimatti et al. (2002) who obtained spectroscopic 
classifications for 29 EROs with K < 19.2 and (R-K) > 5 from 
a full sample of 45. The median redshift of their K < 19.2 sam- 
ple is z = 1 . 1 ± 0.2, with galaxies spread across z = 0.7-1 .4 (see 
Daddi et al. 2001) and they find that their sample is split equally 
into passive, evolved and active, dusty EROs: 31-64% versus 
67-33%, where the ranges reflect the uncertainty due to the in- 
completeness in their survey. The K20 survey has significantly 
increased our knowledge of the ERO population, however, it is 
reliant on identifying emission lines in the restframe UV spec- 
tra of possibly dusty galaxies, which results in an incomplete 
and perhaps biased view of the mix of galaxies within the pop- 
ulation. Moreover, the magnitude limit achieved for the high 
completeness sample, K < 19.2, means that it can't be used to 
track the variation in the properties of the population across the 
break in the counts of EROs. 

Here we apply a new technique to investigate the mix of the 
ERO population at faint magnitudes: using very deep radio ob- 
servations at 1 .4 GHz which are both sensitive enough to iden- 
tify strongly star forming galaxies out to high redshifts, z Z 1-2, 
and yet insensitive to the effects of dust obscuration (Mohan et 
al. 2002; Chapman et al. 2002a). Radio surveys to sub-mJy flux 
limits have shown the usefulness of this technique for identify- 
ing star-forming galaxies in the distant Universe - morpholog- 
ical classifications from Hubble Space Telescope (HST) imag- 



ing show many are blue disk galaxies (Windhorst et al. 1994; 
Richards et al. 1998; Richards 2000), with optical spectroscopy 
confirming that these are apparently normal star forming galax- 
ies at intermediate redshifts (Benn et al. 1993; Mobasher et al. 
1999; Roche, Lowenthal & Koo 2002), with a small proportion 
of AGN (20% based on their radio morphologies in very high- 
resolution maps from combined VLA and Merlin observations, 
Muxlow et al. 1999, with most of these having blue colors). 
However, the spectral properties of the galaxies in these sur- 
veys also show signatures of the presence of dust, which may 
be cloaking the strength of the activity in these systems and 
making them appear more mundane than they should (Hammer 
et al. 1995; Poggianti & Wu 2000; Smail et al. 1999b). 

More recently, samples of faint radio sources have begun to 
be used to attempt to identify obscured star-forming galaxies at 
even higher redshifts. By focusing on those radio sources with 
faint or undetectable counterparts in the optical these surveys 
attempt to isolate the most distant star forming galaxies in the 
radio samples (Richards et al. 1999; Barger et al. 2000; Chap- 
man et al. 2001, 2002b). Sensitive submillimeter observations 
of these optically-faint radio sources confirm that a large frac- 
tion of them are luminous, dusty systems at redshifts of z ~ 1- 
3 - although, as with the whole submillimeter galaxy popula- 
tion, there is still considerable uncertainty over their exact red- 
shifts (Chapman et al. 2002b). The relationship between these 
highly obscured, but very luminous galaxies and the perhaps 
less-obscured systems selected through Lyman-a emission and 
the Lyman-break technique is a crucial issue for a complete un- 
derstanding of the formation and evolution of galaxies (Smail 
et al. 2002, 2003). The project described here seeks to trace 
the population of the luminous, dusty galaxies detected in the 
submillimeter at z k, 2 to lower redshifts and lower luminosi- 
ties. These samples should be more amenable to detailed study 
and so provide a clearer connection between the evolution of 
the obscured and unobscured galaxies over the lifetime of the 
Universe (Chapman et al. 2002b). 

In this paper we use an extremely deep VLA radio map to 
investigate the radio emission from a sample of ERO galaxies 
at a sensitivity limit sufficient to detect strongly star forming 
galaxies at z > 1-2. We also present extensive optical/near- 
infrared imaging of this ERO sample. We exploit this deep 
RIzJHK photometry to analyse the spectral energy distribu- 
tions (SEDs) of the EROs and investigate their properties in 
detail. We present our observations, their reduction and cata- 
loging in the next section, describe our analysis and results in 
§3, and discuss these in §4 before giving our conclusions in §5. 
Throughout we assume a cosmological model with q = 0.5 and 
H„ = 50 km s" 1 Mpc" 1 . If we instead adopted the currently fash- 
ionable H„ = 70kms _1 Mpc -1 , Q\ = 0.7, £l m = 0.3 cosmology, 
then for a source at z = 1 all the linear sizes we derive would 
be 6% smaller, the luminosities would be 12% fainter and the 
volume densities would be 12% lower. 

2. OBSERVATIONS, REDUCTION AND CATALOGING 

The field used for our analysis contains the z = 0.41 cluster 
A 851 (CI 0939+4713). This well-studied cluster is relatively 
rich, but very irregular in appearance (Dressier et al. 1994, 
1997, 2002; lye et al. 2001; Schindler et al. 1998). As we 
discuss below, although cluster members may contribute sig- 
nificantly to the optically-bright radio population in this field, 
the extremely red radio sources discussed in this work all prob- 
ably lie at z > 1 and hence are unassociated with the foreground 
cluster. 
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The other influence which the cluster may have is as a gravi- 
tational lens. However, A 851 is not a particularly concentrated 
cluster (Seitz et al. 1996; Trager et al. 1997; lye et al. 2000) 
and the wide field coverage of our sample, compared with the 
cluster's critical radius (only two of the EROs in our sample lie 
within 1 ' of the cluster center, roughly corresponding to 2-3 x 
the critical radius), suggests that lensing by the cluster will have 
a minimal effect on the surface density in our sample (< 5% 
integrated over the whole area of our survey) and is unlikely 
to grossly effect the morphologies of any of the EROs in our 
study. 

2.1. Radio Observations 

The radio image of A 851 comprises a combination of A, 
B, C, and D configuration observations from the National Ra- 
dio Astronomy Observatory's VLA 9 taken at 1.4 GHz between 
1996 and 2000. For each observation, "4" mode was used for 
the correlator with 7 x 3.125 MHz spectral channels for each 
of two IFs and two polarizations. This mode was employed to 
allow a map to be made of the entire primary beam at the high- 
est spatial resolution possible with the VLA (~ 1.5" FWHM) 
and with the best compromise between radial smearing of the 
image due to the finite bandwidth and sensitivity. However, as 
is noted below, some penalty does result from this choice. The 
total integration time of the combined dataset is ~ 100 hours on 
source. 

The dataset was imaged using the MPS program IMAGR. The 
74 facets were used including 37 within the primary beam and 
37 outside the primary beam, each of the latter centered on a 
bright confusing source. Each of the datasets from the four 
configurations were self-calibrated using intermediate images 
as input to the AIPS program CALIB. A more complete descrip- 
tion of the reduction process is included in Owen et al. (2002). 

The analysis presented here focuses on the central facet from 
the full map, a 8.53' x 8.53' field in the center of the primary 
beam: 09 42 48.614 +46 59 59.72 (J2000). The final map 
reaches a 1-a noise level of 3.5 /zJybeam" 1 in the central part 
of this field - making this the deepest 1.4-GHz map ever made. 
However, in the extreme corners of the area being analyzed the 
actual peak flux density is reduced by about 20% due to the pri- 
mary beam attenuation and radial bandwidth smearing caused 
by the finite bandwidth of each spectral channel. We therefore 
adopt a conservative 3-a limit of 12.6 ^Jy in the corrected peak 
flux density over the entire area being studied. We can probe 
down to 3-a significance limit and still retain a very low level 
of false detections as we are searching for radio emission at the 
known positions of a relatively small number of EROs, com- 
pared to the number of independent beams in our map (10 5 ). 
The synthesized beam of our final map is circular and has a 
FWHM of 1.6", this allows us to resolve radio sources with 
sizes > 1" and we discuss the sizes of radio emission from the 
EROs in our survey in the following sections. 

1 

2.2. Near-infrared Observations 

Our primary near-infrared imaging dataset comes from wide- 
field images in JHK taken with the KPNO 2.1-m telescope, 
these cover a 15.4' x 15.7' field of view to K ~ 21. These im- 
ages have relatively coarse sampling and so we have also ob- 



tained higher resolution J and K images covering a smaller, 
7.5' x 7.5', field of view from the 4.2-m William Herschel 
Telescope (WHT) 10 and the 3.8-m UK Infrared Telescope 
(UKIRT) 11 respectively. The WHT and UKIRT observations 
cover the bulk of the region included in the analysis here and 
they thus provide a useful confirmation of the identifications 
and photometry of EROs in this field. Standard reduction pro- 
cedures were applied to all datasets and they were calibrated 
using UKIRT Faint Standards (Hawarden et al. 2001). 

2.2.1. KPNO 2.1 m JHK-band mosaics 

The SQIID camera was used at the KPNO 2.1-m telescope on 
the nights of 2001 January 4-11 to image the field simultane- 
ously in JHK. The pixel scale of the camera is 0.68 arcsec and 
the useful field of view of is <~ 5.2' x 5.2', with nine pointings 
mosaiced together to cover a region of 15.4' x 15.7'. The total 
exposure time was 75 ks, spread equally over the nine pointings 
to give a per-pixel exposure time of 8.3 ks. 

The images were reduced in a standard manner using IRAF 
scripts. Flat fields were constructed from the dithered expo- 
sures at each pointing to produce a flat field centered in time 
on each exposure. Distortion corrections were measured from 
a comparison of bright objects in the near-infrared images with 
existing, deep /?-band images and these were applied to align 
the images in the three passbands to a common coordinate grid, 
at a level of 0. 1 pixel rms or better. 

These observations were obtained in a mix of photomet- 
ric and non-photometric conditions, with seeing of 1.5-1.8" 
FWHM. To calibrate the photometry in the three bands we ex- 
ploit the precise JHK photometry of this field published by 
Stanford, Eisenhardt & Dickinson (1995). By comparing the 
photometry for a large number of bright galaxies in the field 
we estimate a typical zero-point error of < 0.03 mags on our 
photometric scale. The total magnitude corresponding to the 
5-a detection limit for a point source is K = 20.5. While the 
3-a limits for 3" aperture photometry in the J- and //-bands are 
7 = 22.7 and// = 21.9. 

2.2.2. UKIRT K-band mosaic 

On the nights of 1999 February 28, March 01-02 and March 
15-16 we completed a ZT-band mosaic of 5 x 5 pointings using 
the UFTI imager on UKIRT (pixel scale 0.0908" pixel" 1 ). This 
mosaic covers the central 7.5' x 7.5' region of the cluster de- 
fined by the WFPC2 field (see §2.3) and hence covers around 
55% of the radio sample analysed here. Each pointing has a to- 
tal integration time of 2.4 ks and reaches a limit of K ~ 20.0 for 
point sources. The median seeing for these observations was 
0.5" and several of the mosaic tiles have 0.35" seeing, while 
periods of sub-0.25" seeing were experienced. These frames 
allow us to confirm the reality of EROs in the field and investi- 
gate the morphologies of the brighter examples. 

2.2.3. WHT J -band mosaic 

On 2000 December 11-12 we used the INGRID near- 
infrared imager (Packham et al. 2002) on the WHT to obtain 
a /-band image of the cluster core. This comprised a 2 x 2 
mosaic of 2.4-ks exposures and covers the same region as the 
UKIRT mosaic, reaching J = 22.7 for a point-source in 0.70" 
seeing with 0.242" sampling. 



9 NRAO is operated by Associated Universities Inc., under a cooperative agreement with the National Science Foundation. 
10 The WHT is operated by the Isaac Newton Group on behalf of Particle Physics & Astronomy Research Council (PPARC). 
1 1 UKIRT is operated by the Joint Astronomy Centre on behalf of PPARC. 
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Fig. 1. — The (R-K)-K color-magnitude diagram for 1186 sources brighter than K = 20.5 within a 8.53' x 8.53' field centered on A 851 - these 
are shown as small points and lower limits. Overlayed on this are the radio-detected sources brighter than the 4-a limit of our radio map, 16.8 /ijy. 
These radio sources go fainter than the nominal K = 20.5 catalog limit (shown by the dotted line). The solid line shows the (R-K) = 5.3 boundary 
used to define our ERO sample. For the ERO sample we identify all sources detected above a slightly deeper radio limit, 12.7 /ijy or 3a, again 
indicating their radio fluxes by the size of the symbols used. The dashed line denotes the reddest color expected for unobscured galaxies at the 
cluster redshift; galaxies redder than this limit are likely to be beyond the cluster. The population of radio-emitting cluster members is clearly 
visible below this boundary (Morrison 1999; Morrison et al. 2002). 



2.3. Optical Observations 

The optical imaging to complement our deep near-infrared 
datasets was obtained from three facilities. We have two wide- 
field imaging datasets of this region: the first comes from 
Suprime-Cam on Subaru (courtesy of Dr T. Kodama), and con- 
sists of very deep RI imaging of a 27' x 27' field centered on the 
cluster core. We also have obtained similarly wide-field KPNO 
4-m image in the z-band to bridge the spectral coverage of our 
optical and near-infrared datasets. Finally, an HSTAVFPC2 mo- 
saic of A 85 1 has been observed in the F702W filter, this covers 
a 7.5' x 7.5' region around the cluster core, within the area stud- 
ied here, and provides very detailed morphological information 
on the optically-bright galaxies in this area (Dressier et al. 2002; 
Morrison et al. 2002). However, the relatively shallow surface 
brightness limit of the HST imaging, /ir ~ 25.5 mag arcsec -2 , 
means that these data are not useful for gauging the morpholo- 
gies of all but the most compact of the 15 radio-detected and 27 
radio-undetected EROs in this region (see also Dey et al. 1999; 
Smith et al. 2002a). For this reason we will not use the HST 
imaging in our analysis. 

2.3.1. Subaru RI imaging 

Kodama et al. (2001) obtained deep, RI imaging of A 85 1 us- 
ing Suprime-Cam on the 8.3-m Subaru Telescope on the nights 
of 2001 January 21-22. Suprime-Cam has a 27' x 27' field 



of view, imaged with eight 4kx2k CCDs at a scale of 0.20" 
pixel" 1 . The conditions were photometric during the observa- 
tions, with seeing of 1 .0" and 0.7" in R and / respectively. Total 
exposure times of 4.0 ks and 1.3 ks were obtained resulting in 

3- <7 limiting magnitudes for point sources in 3" diameter aper- 
tures of R = 26.7 and / = 25.6. The reduction and calibration of 
these data is described by Kodama et al. (2001). 

2.3.2. KPNO 4m z imaging 

The z-band imaging of A 85 1 was obtained on the night of 
2001 January 25-26 with the MOSAIC camera on the KPNO 

4- m. These data cover a 36' x 36' field with 0.258" pixel -1 
sampling. The processed image combined five individual expo- 
sures, shifted so the combination eliminated areas lost to gaps 
between individual CCD chips. After applying chip-by-chip 
sky flats, each exposure was projected onto an astrometric pixel 
grid and coadded with appropriate masking for gaps and bad 
pixels. The data were taken in photometric conditions and the 
final 7.2 ks exposure has seeing of 1 .0" FWHM and a 3-cr limit- 
ing magnitude of z = 24.1 in our standard 3" diameter aperture. 

2.4. Photometric Catalog 

To catalog near-infrared sources in the field of A 85 1 we ran 
SExtractor (Bertin & Arnouts 1996) on the K-bwd image of the 
field - first convolving the image with a 2"-diameter top-hat 
kernel and then searching for sources with areas greater than 
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the seeing-disk above an isophote defined as /x^ = 22.8 mag. 
arcsec" 2 . 

As our measure of the total magnitudes of the sources in the 
image we adopt the scaled aperture magnitude, BEST_MAG, 
calculated by SExtractor. This magnitude is measured within 
an aperture which is 2.5 x the first-moment radius of the source, 
except in crowded environments when it is replaced with a cor- 
rected, isophotal magnitude, see Bertin & Arnouts (1996) for 
details. The 5-cr detection limit for point-sources in our cata- 
log corresponds to a total magnitude of K = 20.5 and we detect 
1186 sources in the area spanned by our current VLA catalog 
(Fig. 1). This point-source detection limit should be representa- 
tive of the completeness of our catalog for those compact galax- 
ies with half-light radii less than our 2"-diameter detection filter 
(this includes the majority of the faint, near-infrared field pop- 
ulation, Chen et al. 2002). However, our catalog selection will 
tend to miss galaxies with very extended disks even though their 
integrated magnitudes are above this point-source limit and this 
may introduce a morphological bias into our sample selection 
(Martini 2001; Chen et al. 2002; Firth et al. 2002; Snigula et al. 
2002). However, this is not a significant problem for our anal- 
ysis as the primary goal is to determine the proportion of dusty 
and evolved EROs in published surveys and the relatively good 
agreement between our measured ERO surface density and val- 
ues published by other workers (see Fig. 2) suggests that our 
catalog suffers no stronger bias than is typical for such surveys. 

To measure colors for galaxies in our RIzJHK imaging 
we use fixed, 3"-diameter aperture photometry from seeing- 
matched frames. The optical and near-infrared images were 
aligned to the coordinate system defined by the KPNO TT-band 
image with a typical tolerance of < 0.1". The seeing on these 
frames was then matched to that of the poorest-resolution im- 
ages (the J /H /K-band) by Gaussian convolution. 

3. ANALYSIS AND RESULTS 

Our analysis begins by constructing a SExtractor-based K- 
selected catalog of EROs in the field. We then align the radio 
and optical/near-infrared images to an rms accuracy of 0.32" 
using the 90 brightest radio and TT-band sources in the joint 
field. This allows us to search for matches for the EROs within 
the radio map. We note that our radio observations identify all 
three cataloged submm galaxies in this field (Smail et al. 1999a, 
2002; Cowie et al. 2002) and we discuss these further in Ledlow 
et al. (2002). 

3.1. The ERO samples 

For this analysis we start from our Zf-selected galaxy catalog 
and identify galaxies brighter than K = 20.5 (5a) which have 3- 
" diameter aperture colors redder than the limit of (R-K) = 5.3 
which we adopt as the definition of an ERO (see Pozzetti & 
Mannucci 2000). Based on this we identify 68 EROs in our 
field - equivalent to a surface density of (1 .04 ± 0.12) arcmin" 2 
at K = 20.5 (the error-bar simply reflects Poisson statistics and 
we have applied an incompleteness correction of 20% to the 
number density of EROs in the faintest bin, as well as a 6% 
correction to the area of the survey to account for the sky area 
covered by bright stars and galaxies). 

The catalog of EROs in our field is given in Table 1 . This lists 
the ERO ID (a 5-digit ID denotes a radio source), the J2000 po- 
sition of the ERO (tied to the VLA reference frame with an rms 
precision of 0.32"), the total TT-band magnitude (measured us- 
ing the SExtractor BEST_MAG), aperture (R-K) color (or 3-cr 
limit), 1.4-GHz radio flux (or 3-cr limit), 3-" diameter aperture 



magnitudes from our seeing-matched RIzJHK imaging, best- 
fitting SED (either dusty, "D"; evolved, "E"; or badly-fit, "B", 
see §3.5) and a visual morphology (from our high resolution K- 
and /?-band imaging if available, see §3.4) on a simple system: 
"C", compact/regular; "I", disky, disturbed or amorphous; "F", 
too faint to classify. 

1 i 1 i 1 i 1 

-•Whole sample, (H-K)>5.3 
▲ Radio detected, S, 4 >12.6/ijy 




K 

Fig. 2. — The cumulative number density of EROs in our sample 
compared to published counts from Daddi et al. (2000), Thompson 
et al. (1999) and Smith et al. (2002a). We show both the full ERO 
catalog, selected with (R-K) > 5.3 and also the radio-detected sam- 
ple. The error-bars represent simple Poisson uncertainties and are 
likely therefore to be underestimates of the true field-to-field varia- 
tions (Daddi et al. 2001). We have included a 6% downward cor- 
rection to the estimated area of our survey due to coverage by bright 
stars and galaxies in our field. The dashed line shows the fit to the 
radio-detected ERO population with a normalization increased to re- 
produce the 60% dusty ERO fraction estimated for our whole sample 
to K = 20.5. The dotted line is a Gaussian component with a mean 
magnitude of K = 19. 1 and a FWHM of 1.2 mags which is normalized 
so as to reproduce the cumulative counts of EROs in our field. The 
solid line shows the behavior of a toy model constructed by combining 
the passive and dusty populations indicated by the dashed and dotted 
lines. The dot-dashed line shows a second toy model which combines 
a series of Gaussian luminosity functions for the ERO population with 
a mean of L* K and a dispersion of ok = 1 mag, which are coadded so 
as to produce a Gaussian redshift distribution with <z>~ 1.1 and 
a z = 0.2 (see Firth et al. 2002). 

As is shown in Fig. 2, the surface density of EROs we detect 
is in reasonable agreement with that observed in other ERO sur- 
veys (e.g. Thompson et al. 1999; Daddi et al. 2000; Smith et al. 
2002a) given the strong field-to-field variation seen in the sur- 
face density of EROs (Daddi et al. 2001; McCarthy et al. 2001; 
Roche et al. 2002; Chapman, McCarthy & Persson 2000). This 
supports our claim that the foreground cluster makes no con- 
tribution to this population, either from members or by grav- 
itational lensing. For ease of comparison with ERO samples 
selected in other filter combinations, in particular (I-K), we 
note that while our field contains 68 EROs with (R-K) > 5.3, 
we have 99 objects with (I-K) > 4 suggesting that this is a 
less stringent definition of an unusually red source. We also 
note that 63 (93%) of our (R-K) > 5.3 EROs have (I-K) > 4 
(including all of the radio-detected EROs from §3.2). 
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We re-identify three previously cataloged EROs brighter 
than K = 20.5 in our field, including the submm ERO 
SMMJ09429+4658 from Smail et al. (1999a), which is also 
a faint radio source, and two other radio-undetected EROs 
- the edge-on disk galaxy #333 from the NICMOS study of 
Smail et al. (1999b) and CI 0939+47 13 A from Persson et al. 
(1993). A fourth cataloged ERO in this field - the passive 
z = 1.58 ERO CI 0939+4713 B from Soifer et al. (1999) - has 
(R-K) = 5.1 ±0.1 based on our photometry and hence does not 
make it into our ERO sample. 

We can also use our UKIRT and WHT near-infrared imag- 
ing to test the reliability of our ERO sample. There are 30 
EROs from the KPNO sample which fall within the UKIRT K- 
band and WHT /-band mosaics (1 1 of these are radio-detected). 
We identify 27 (90%) of these on the UKIRT mosaic. Two of 
the missing EROs are low-surface brightness, extended systems 
which are below the detection limit of the UKIRT mosaic when 
using a large diameter photometric aperture. This highlights 
the surface brightness, and hence morphological, biases which 
could affect the mix of morphological classes identified in ERO 
surveys (Chen et al. 2002). In this regard the good sensitivity, 
but relatively coarse, resolution of our KPNO /T-band imaging 
acts in our favor by concentrating a large fraction of the light 
from the galaxy in a relatively small number of pixels (50 x 
fewer than in the UKIRT image). 

The third "ERO" missing from the UKIRT mosaic is the 
source #2549 cataloged in the KPNO 2.1-m imaging (Table 1). 
This source is clearly detected in the simultaneously obtained 
JHK imaging from 2001 January 4-1 1 as a bright K = 19.0 un- 
resolved source (< 1.5" FWHM) with unusual colors (J-K) = 
0.22 ±0.07, (H-K) = -1.09 ±0.07: indicating that the SED 
peaks in the fl-band (Fig. 7). Yet it was apparently fainter than 
K > 20.0 during February 1999 (UKIRT) and J > 22.7 in De- 
cember 2000 (WHT). In addition, we can state the counterpart 
to this source was also fainter than our 3-a limits of R > 26.7 
and / > 25.6 on 2001 January 21 and z > 24.1 on 2001 January 
25. The strongest constraint on the variability of the source 
comes from our /-band observations which would indicate that 
the source brightened by 3.5 magnitudes in a 1 month period. 
While our subsequent optical images suggest it was either ex- 
tremely red, (/-/)> 6.5 or it had faded by > 4 magnitudes 
within 3 weeks of the 2.1-m run. There are two possible ex- 
planations for the properties of this source: it is either moving 
(e.g. an asteroid or comet) or highly variable (e.g. a supernova, 
gamma-ray burst or active galactic nucleus). We can easily dis- 
count at least one of the variability options - if the source was 
a supernova then it should still be visible in our optical imag- 
ing taken 3 weeks after the 2.1-m run. Equally, by subdividing 
our 2.1-m imaging into temporally independent sections and 
comparing the position of the object between these, we find 
no measurable movement of the source, < 0.5", or strong vari- 
ability over a 6 hr timespan. If the object was moving at the 
maximum rate allowed by this limit then it should be still visi- 
ble within our z-band frame taken ~ 3 weeks later. Its absence 
therefore rules it out as a solar-system object. We conclude that 
#2549 is most likely a highly variable source, although unlike 
the highly-variable AGN found by Gal- Yam et al. (2002) there 
is no optical or infrared host galaxy underlying #2549. We also 
note that #2549 lies 30" from the variable X-ray source in this 
field cataloged by Schindler et al. (1998) and so is unlikely to 
be the same object. 

Based on these comparisons we conclude that our 
photometrically-selected ERO sample is reliable. 



3.2. The Radio-detected EROs 

We now search for radio counterparts to the EROs in our 
field. We adopt a maximum matching radius of 2" (3 pix- 
els on the near-infrared images) and search for radio emission 
above the 3-(T limit of 12.6/iJy centered on the Zf-band cen- 
troid of each ERO. Of the 68 photometrically-selected EROs, 
21 are detected in our deep VLA observations - giving a surface 
density of (0.31 ±0.07) arcmin" 2 for EROs with radio fluxes 
above 12.6 yuJy and brighter than K = 20.5. The median TT-band 
magnitude of these 21 EROs is K = 19.66 ± 0.22 (Table 2), al- 
most identical to the median brightness of the radio-undetected 
EROs: K = 19.65 ± 0.09. The brightest ERO in the radio has 
a flux of only 200/iJy, we can therefore limit the proportion of 
radio luminous AGN in the general ERO population to < 2% 
(c.f. Willott et al. 2001). 

The radio-detected EROs comprise 31 ±7% of the whole 
ERO population down to K = 20.5 (see also Yan 2001 ; Mohan et 
al. 2002), where the error-bar is based on Poisson statistics and 
hence probably underestimates the field to field variance in this 
fraction. Taking the expected blank-field radio counts sources 
down to 12.6 /iJy (<~ 250, Richards 1999) the ERO counterparts 
with K < 20.5 represent ~ 8% of the total radio population. 

We plot the cumulative number density of the radio-detected 
EROs in Fig. 2. We find, at most, marginal evidence for a 
difference in the count slope between the whole ERO popu- 
lation (and given our small sample also no strong evidence 
for a break in the count slope with magnitude) and the radio- 
detected and undetected samples: parameterizing the counts as 
log 10 (A^ERo) oc aK we estimate a = 0.59 ± 0.05 for the whole 
catalog, with a = 0.5 1 ± 0.08 and a = 0.64 ± 0.07 for the radio- 
detected and undetected samples. 

The radio flux distribution of the radio-detected EROs is 
characterized by a strong rise at faint fluxes - the cumulative 
counts being well-fit by a S] 4 19±om power-law. Hence the 
number of EROs detected at 1 .4 GHz rapidly increases at fainter 
flux limits: doubling as the flux limit is reduced from ~ 40/iJy 
(typical of previous surveys, e.g. Richards et al. 1998) to the 
12.6-ytiJy limit of the current work. At a flux level of around 
50 /iJy the ERO counts appears to increase more rapidly (per- 
haps indicating the appearance of a new population of very red 
and very faint radio sources), but then flattens out at faint fluxes. 
We believe that this flattening is probably a result of the depth of 
our available TT-band imaging which is restricting the identifica- 
tion of the faintest and typically reddest sources (see Fig. 3). We 
see fairly constant /f-band to radio flux ratios as a function of 
radio flux (see Fig. 3), this is in contrast to the strong decrease 
in this ratio with increasing radio flux above 1 mJy (Waddington 
et al. 2000), suggesting a different physical process (star forma- 
tion, rather than AGN activity) may be responsible for the faint 
emission from the sources in our catalog. 

The distribution of EROs in Fig. 3 shows that these galaxies 
have typically lower near-infrared to radio flux ratios than aver- 
age radio source in the field (many of which are cluster mem- 
bers, Morrison et al. 2002). The similarity of the K-corrections 
in these two wavebands over the redshift range we cover (see 
§3.5) suggests that this is an intrinsic feature of these galaxies. 
Assuming that the local far-infrared to radio correlation holds 
for these galaxies, then we conclude that the radio-detected 
EROs have far-infrared/optical ratios of L FIR /L V <~ 20-200 (see 
§4.2). Locally such high Lfju/Ly ratios are usually only seen 
in the advanced stages of mergers of massive galaxies (e.g. 
Mirabel & Sanders 1989). In comparison, using the median 
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TT-band magnitude for the radio-undetected EROs we obtain 
a limit of Sk/S\a > 0.6, equivalent to L FIR /L V £ 25 (for an 
Ly star-forming galaxy at z > 1 the equivalent FIR luminosity 
would translate into a limit on the SFR of > 0.1 M Q stars of 
< lOOMoyr" 1 ). 




-10 100 1000 



s, 4 C"Jy) 

Fig. 3. — The S K /Siaghz-Siaghz plane illustrating the variation in 
near-IR to radio flux with radio flux for all the radio-detected sources 
in our field. The dashed line shows the selection boundary for sources 
brighter than K = 20.5 and 12.6 pjy and the bold symbols identify 
those radio sources which have ERO counterparts. We also illus- 
trate the tracks for SEDs with far-infrared to optical flux ratios of 
Lfir/Lv ~ 5 and 50 (using the far-infrared/radio correlation for lo- 
cal galaxies, Condon 1992). We mark on these tracks the equivalent 
redshift limits (we mark 8z = 0.5 increments on each track) assum- 
ing a model where the source is an Ly(M v = -21.8) galaxy with an 
SED similar to a present-day Sbc galaxy and radio luminosities of 
U a = 2.3 x 10 23 (upper track) and 2.3 x 10 24 WHz 1 (lower) with a 
radio spectrum of the form 5 oc v^' 1 . These radio powers are equiv- 
alent to SFRs of ~ 20 and ~ 200 M yr" 1 for stars above 5 M Q , re- 
spectively, or equivalently 10 2 and 10 3 Mq yr~' assuming a Salpeter 
IMF which extends down to 0.1 Mq (see §4.2). 

Finally, we searched the archival 49.4-ks XMM- 
Newton/EPIC X-ray image of this field to identify any ERO 
which may be associated with a luminous X-ray sources 
(Alexander et al. 2002; Stevens et al. 2002; Page et al. 
2002). We find no such X-ray bright EROs within our 72.7 
sq. arcmin search region and place typical 3-a limits on 
the unabsorbed X-ray emission from individual sources of 
< 4 x 10" 16 ergs s" 1 cm" 2 in the 2-10 KeV band. 

3.3. ERO companions 

To search for differences in the local environment of 
the radio-detected and undetected ERO populations we look 
at the typical nearest-neighbor distance in our K < 20.5 
photometrically-selected catalog (Georgakakis et al. 1999; Gal- 
limore & Keel 1993). We note that the modest spatial resolu- 
tion of the KPNO /T-band images means this test is relatively 
insensitive to the presence of close companions, < 1-2" seper- 
ation. We find that the median nearest-neighbor distance is 
10.9"±1.6" for radio-detected ERO, compared to 15.2"±0.8" 
for those EROs which are not detected in our radio map. A 



Kolmogorov-Smirnov test shows only a 10% chance that the 
two nearest-neighbor distributions are drawn from the same 
parent population. This suggests a higher frequency of neigh- 
bors for those EROs detected in our radio map, with sepera- 
tions roughly corresponding to ~ lOOkpc at z > 1. However, 
these seperations are probably too large for the companion to 
be responsible for dynamically triggering the activity in the ra- 
dio source. To investigate the sub-10-kpc scales, where the sig- 
natures of on-going mergers and strong dynamical interactions 
will be visible, we need to look at the detailed morphologies of 
the EROs (§3.4). 

To further test the association of these neighboring galax- 
ies with the EROs we analyse their colors. We find that the 
radio-detected sample has a median nearest-neighbor color of 
(R-K) = 3.48 ± 0.20, while the closest K < 20.5 galaxy to 
the radio-undetected EROs have (R-K) = 3.96 ±0.31 - in 
both cases the nearest neighbors have a median magnitude of 
K ~ 19.2 ± 0.2. Both of these colors are marginally redder than 
the average for a magnitude-matched subsample from the field: 
(R-K) = 3.26 ± 0.05, suggestion that the nearest neighbors to 
the EROs may lie at a higher mean redshift, are more evolved 
or more obscured than the typical field galaxy. Both samples 
also have comparable fractions of EROs as nearest neighbors 
- 14% of the radio undetected sample have a nearest neighbor 
which is also an ERO (with (R-K) > 5.3) compared to 12% for 
the radio-detected EROs. Given that EROs make up only 6% of 
the whole /f-band sample in this field it is clear that the EROs 
as a population must be clustered (Daddi et al. 2001; McCarthy 
et al. 2001 ; Chapman, McCarthy & Persson 2000). 

3.4. ERO morphologies 

We can also exploit our high-resolution UKIRT /T-band and 
Subaru /?-band imaging to quantify the morphological mix 
within the two samples of EROs (Fig. 4). The UKIRT imaging 
provides restframe optical morphologies for 19 of the radio- 
undetected EROs (including the transient source #2549, §3.1, 
and the star #2418, §3.5) and 11 radio-detected EROs. The 
properties of these subsamples are consistent with the full sam- 
ples (Table 2) with median TT-band magnitudes and (J-K) col- 
ors of 19.61 ±0.16, 2.09±0.12 and 19.90±0.24, 2.23±0.12 
respectively. In contrast the Subaru /?-band image covers the 
full ERO sample at slightly lower resolution, but only gives us 
a restframe UV view of the morphologies of these galaxies - 
which is very sensitive to the presence of both dust and unob- 
scured star formation. 

We visually classify the EROs in both passbands using a 
simple scheme designed to crudely distinguish between regu- 
lar, centrally-concentrated galaxies, whose extreme colors may 
arise from a old stellar population, and more morphologi- 
cally complex, disturbed or asymmetric galaxies, where their 
very red colors could indicate the presence of dust obscura- 
tion. Following Smith et al. (2002a) we define three broad 
classes: compact/relaxed systems ("C"); obvious disks, asym- 
metric or disturbed/amorphous morphologies ("I") or too faint 
to reliably classify ("F"), see Table 1. Based on the near- 
infrared classifications we estimate that the morphologies of 
the radio undetected ERO sample are: 5 compact, 5 disks 
or disturbed/amorphous and 7 which are too faint to classify 
(29%:29%:42%), excluding #2549 and #2418. For the radio- 
detected sample the equivalent numbers are: 0, 9 and 2 respec- 
tively (0%:82%:18%). In both cases the uncertainties in the 
numbers in each class are probably 1-2 objects. We confirm 
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Fig. 4. — Example of the R- and A-band morphologies of a random selection of radio-detected (top two rows) and undetected (bottom two rows) 
EROs which fall within the UKIRT mosaic. Each field is shown in two panels, with the Subaru i?-band images above the UKIRT A"-band in both 
cases. We label the panels with the catalog number, AT-band morphology and radio flux (if detected). There is a clear tendency for radio-detected 
EROs to show more complex morphologies, while radio undetected EROs are more regular and compact. We show two sources in the panel 
containing the compact ERO #2176 and an example of an amorphous, radio-undetected ERO: #2128. Note that #1931 is the edge-on disk galaxy 
identified in the NICMOS image from Smail et al. (1999b) and #01208 is the ERO counterpart to SMM J09429+4658 (Smail et al. 1999a). Each 
panel is 12" x 12" square and has North top and East to the left, the typical seeing is 0.5" in the A"-band and 1.0" in the A-band. 
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that these two distributions are inconsistent with being drawn 
from the same parent population at the 10 -4 level. 

The morphological mix of the whole sample is: 
18%:50%:32%, identical to that measured for a K < 21 .5 lensed 
sample by Smith et al. (2002a). We find a similar distribution 
when we use the restframe UV classifications from our /?-band 
images: with the mix for the full sample being 1 1%:52%:37%, 
and the radio-detected EROs having a slightly higher incidence 
of disturbed/amorphous morphologies and a lower incidence of 
compact images than the radio-undetected subsample. Thus our 
sample appears to exhibit a lower fraction of compact, regular 
EROs compared to that presented by Moriondo et al. (2000), 
where 80% of their EROs were classified as compact/elliptical- 
like from an analysis of a somewhat brighter (median K = 19.3) 
archival HST sample. This difference could either reflect a mis- 
match in the classification schemes; a tendency for the archival 
HST fields to target high-density regions, which may contain 
a higher fraction of bulge-dominated EROs; or a real change 
in the morphological mix of the ERO population with apparent 
magnitude (see §4.2). 

Turning to the internal differences within the subsamples of 
EROs in our survey, although the samples are small and the 
classifications are necessarily crude and somewhat uncertain, it 
appears that the radio-detected EROs show a higher incidence 
of disky and disturbed/amorphous morphologies compared to 
the radio-undetected sample, which in contrast has a higher 
fraction of compact sources (no examples of which are seen 
in the radio-detected sample). We illustrate a few examples 
of the various class of ERO morphologies in Fig. 4. We also 
briefly note the median (J-K) colors of the different morpho- 
logical classes (both radio-detected and undetected): compact, 
(J-K) = 2.02 ± 0.36; disky/amorphous, (J-K) = 2.32 ± 0.15; 
faint, (J-K) = 2.09 ±0.11. Thus there is a hint that the com- 
pact population is bluer than the extended (disky, disturbed and 
amorphous) subsample, in agreement with the findings of Man- 
nucci et al. (2002), but the large scatter in the colors of these 
subsamples suggests that they do not represent homogeneous 
populations. 

Finally, we note that the radio emission from 7 of the radio- 
detected EROs, 33% of the sample, is resolved in our map on 
scales of 0.8-3.4". The median intrinsic size of these sources is 
1 .4 ± 0.4", corresponding to 12 ± 3 kpc at z ~ 1 in our adopted 
cosmology. The remaining 14 radio-detected EROs are unre- 
solved with a typical limit on the size of their radio emission 
of < 0.8" (< 7 kpc at z ~ 1). For the 7 resolved EROs we can 
crudely compare their morphologies in the radio and optical 
passbands by determining the median offset in the position an- 
gles of the emission in these two bands. We find that the radio 
and /?-band emission have very similar scale sizes and that there 
is a median offset of only 1 1 ± 14 degrees between the position 
angles of the major axes in the two bands. This suggests that 
there is a close relationship between the extended radio and op- 
tical emission in these galaxies and we return to this point in 
§4. 

3.5. ERO Colors 

We begin our investigation of the colors of the EROs in this 
field by using the most basic diagnostic diagram proposed for 
this population: the (R-K)-(J-K) plane. Pozzetti & Mannucci 
(2000) proposed that evolved and dusty EROs could be sepa- 
rated using this plane - with the passive, evolved EROs typi- 
cally showing bluer (J-K) colors (see also Smith et al. 2002a; 
Mannucci et al. 2002). We therefore plot the distribution of all 



the well-detected EROs (see below for the definition of well- 
detected) on this plane in Fig. 5. Focusing on the relative distri- 
bution of the radio-detected and undetected EROs we see that 
the radio-detected EROs tend to fall red-ward, or close to, the 
division proposed by Pozzetti & Mannucci (2000). In contrast, 
the radio-undetected EROs show a much broader distribution, 
with roughly equal numbers of galaxies either side of the divi- 
sion (see also Mannucci et al. 2002). Figure 5 also allows us to 
classify #2418 as a probable star (this is confirmed by the form 
of its SED in Fig. 6). 

To investigate the photometric properties of the ERO popu- 
lations in more detail we now exploit our multiband photome- 
try. We begin by simply constructing average colors for each 
sample. The median colors of the radio-detected and radio- 
undetected ERO samples are given in Table 2 (where the er- 
rors are the uncertainty on the median from bootstrap resam- 
pling). We find that the radio-undetected EROs are typically 
bluer than the radio-detected sample in colors long-ward of 
~ 1 /im - consistent with the color distribution of the two sam- 
ples in Fig. 5. We also find that the radio ERO population may 
be slightly more homogeneous than either the whole ERO cat- 
alog or the radio-undetected sample. For example, the radio- 
undetected sample exhibits a broader spread of colors with a 
typical dispersion of 0.44 ±0.03 around the median color (Ta- 
ble 2), compared to 0.37 ± 0.03 for the radio-detected sample 
(this is also illustrated by the broader range of SEDs for the 
radio-undetected sample in Fig. 6). The varied mix of colors 
for the radio-undetected EROs is not surprising as we expect 
that this sample probably contains both passive evolved galax- 
ies and dusty active systems whose radio flux is below the sen- 
sitivity limit of our VLA map. The presence of passive, evolved 
galaxies in this sample is further supported by the more abrupt 
change in the colors of this sample at < 1/im, corresponding 
to the 4000A break at z < 1.5. However, the relatively mod- 
est differences between the median colors of the radio detected 
and undetected samples also suggests that the radio-undetected 
ERO sample contains a reasonable fraction of galaxies with 
dusty SEDs (or that the relative redshift distributions of the two 
samples conspire so that the observed colors are more similar 
than they are in the restframe). 

3.5.1. Photometric classification of EROs 

To further investigate the photometric classifications of in- 
dividual ERO we have undertaken a simple test to determine 
whether the colors of a particular ERO are better described by a 
relatively feature-less, dust-dominated SED or a spectrum with 
a strong break. The simplest way to apply this test is to em- 
ploy the HYPERZ photometric redshift code (Bolzonella et al. 
2000) to compare the ERO colors to two different families of 
SEDs: dusty, young starbursts (or AGN) and almost dust-free 
evolved systems. For the dusty AGN we assume that their opti- 
cal and near-infrared colors are dominated by a co-eval starburst 
(which produces the dust) - although the AGN may contribute 
significantly to the radio emission from the system. By relying 
on just two simple classes of SED we will clearly fail to de- 
scribe any systems with more complex mixes of evolved stars 
and obscured and unobscured young stars (a situation which 
may be common place in many of these systems at high-z, e.g. 
the edge-on bulge-strong disk galaxy ERO #1931, Fig. 4 and 
Table 1). However, by concentrating on those EROs which are 
well-fitted by these simple models we can determine the rela- 
tive proportion of EROs at the extremes of these classes. 
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Table 2 

Median Colors of ERO Subsamples 



Sample 


N 


K 


(R-K) 


(I-K) 


iz-K) 


(J-K) 


(H-K) 


Dispersion 


Note 




Radio detected 
All 


21 


19.66 ±0.22 


5.88±0.13 


4.53 ±0.15 


3.56 ±0.08 


2.05 ±0.09 


1.18 ±0.08 


0.37 ±0.03 






Radio undetected 
All 


45 


19.65 ±0.09 


5.82 ±0.11 


4.67 ±0.14 


3.40 ±0.07 


1.98 ±0.06 


1.05 ±0.10 


0.44 ±0.03 


Excluding #2418 t 


i #2549 


Dusty SED 
Evolved SED 
Badly fit 


9 
11 

8 


19.67±0.10 
19.60±0.11 
19.22 ±0.24 


5.73 ±0.19 
5.59 ±0.17 
5.73 ±0.32 


4.51 ±0.13 
4.89 ±0.20 
4.27 ±0.22 


3.39±0.18 
3.49 ±0.11 
3.60 ±0.37 


2.53 ±0.17 
1.97 ±0.14 
1.81 ±0.12 


1.21 ±0.11 
1.18=1=0.12 
0.83 ±0.04 


0.36 ±0.03 
0.41 ±0.07 
0.49 ±0.15 


Excluding #2418 <! 


i #2549 


Radio detected & 
Not fit 


undetected 
20 19.92±0.12 


5.81 ±0.13 


4.94 ±0.17 


4.48 ±0.45 


2.19±0.19 


1.15 ±0.15 









For this test we employ a single star formation model de- 
scribed by an exponential declining star formation rate with a 
e-folding time (r) of 1 Gyr, a Miller-Scalo IMF and solar metal- 
licity. This model can provide an adequate description of the 
colors of luminous elliptical galaxies at the present-day (Bol- 
zonella et al. 2000) and hence is a reasonable choice given the 
expectation that both the passive and star-forming ERO popula- 
tions may represent the precursors of the elliptical galaxies seen 
in the local Universe. As a test of the sensitivity of our classifi- 
cations to the details of the star formation model we have also 
used an instantaneous burst and a r = 2 Gyr model in place of 
the r = 1 Gyr model and find that the proportions of EROs in 
the different classes varies by less than 10%. 

In our model the initial formation redshift, the dust redden- 
ing and the redshift of the ERO are all free parameters subject 
to the following constraints. For the starburst/AGNs the red- 
dening, based on a Calzetti et al. (2000) dust model, can vary in 
the range A v = 1-6, consistent with the range of reddening esti- 
mated for their spectroscopically-identified star forming EROs 
by Cimatti et al. (2002). In contrast for the passive models it is 
required to be Ay < 0.5 - hence by restricting the range of pos- 
sible reddening we isolate the two families of solutions which 
are capable of fitting the very red colors of the EROs: either un- 
obscured but old stellar populations or highly obscured, actively 
star forming systems. In both cases the redshift of the ERO is 
allowed to vary between z = 0-3, a large enough range to com- 
fortably encompass all of the EROs in Cimatti et al. (2002). The 
upper limit on the redshift removes the possibility of fitting in- 
feasibly high redshifts due to the 4000A break being confused 
with the Lyman break, although we note that at least one ERO 
in our sample has a proposed redshift of z > 3.4 (#01206, Smail 
et al. 2002). 

To obtain useful constraints from our analysis, we first re- 
quire that the ERO is detected in at least 3 passbands at more 
than 3-cr, this reduces the samples to 17 radio-detected EROs 
and 30 radio undetected EROs. We list at the bottom of Ta- 
ble 2 the median /f-band magnitude and colors of those EROs 
which were removed from our analysis by this signal-to-noise 
ratio (SNR) cut. As expected these are fainter on average, 
K= 19.92 ±0.12, than the typical ERO in our sample, which 
has K = 19.63 ± 0.09. This magnitude difference is equivalent 
to that expected for a change in the median distance modulus 
of the population of 5z ~ 0.1. As many of these EROs are un- 
detected in a particular passband we have improved the signal- 
to-noise of their colors by calculating these from the median 
fluxes in the various passbands, including cases where the cal- 
culated flux is negative (Table 2). We find that these colors are 



consistent with those for the full samples, with the largest dis- 
crepancy being the (z—K) color which differs by ~ 2.2a. We 
conclude that the EROs removed by our SNR cut are at most 
slightly redder than the sample we retain, and the main reason 
for their lower SNR is that they are simply fainter on average. 

We fit the RIzJHK photometry (including limits) of the high- 
SNR radio-detected and radio-undetected samples to both fami- 
lies of model SEDs allowing for dm = 0.05 systematic errors on 
the calibrations of the different passbands. We then compare 
the relative \ 2 of the best-fit models in both cases and classify 
the EROs on this basis into either evolved, "E", dusty, "D", or 
badly-fit, "B" (see Table 2). We plot the SEDs of the sources 
used in the analysis in Fig. 6. Encouragingly this test identifies 
#2418 and #2549 as pathological SEDs which are not well fit 
by either galaxy model (x 2 3> 100) and we remove these from 
our sample, reducing it to 28 sources (Table 2). 

Restricting our comparison to those EROs which are reason- 
ably well-fit by either SED, \ 2 < 2.7 (90% confidence limit) we 
find that 12 from the 17 EROs in the high-SNR, radio-detected 
sample are well-fitted, 1 1 of these by the dusty SED (65%). 
This supports the suggestion that the bulk of this population 
represent dusty, active galaxies and gives us good confidence 
that this model SED can be reliably used to identify dusty sys- 
tems in the high-SNR, radio-undetected sample. We find only 
one radio-detected ERO which is better fit by an evolved SED, 
this is #11077, a relatively faint radio source with an irregu- 
lar ERO counterpart. Examining the SED of this source we 
see that it is unique amongst the radio-detected EROs in hav- 
ing a blue (R — I) color, which cannot be well-fitted by a dusty 
SED (Fig. 6). Comparing the morphology of the ERO in the R- 
and TT-bands we find that the galaxy comprises two regions - 
a blue area to the east of the radio source and a very red area 
to the west - this complex morphology suggests that the simple 
model we use to fit may not be applicable to this galaxy. 

The five high-SNR, radio-detected EROs which are not well- 
fit by the dusty SED (#01045, #01105, #01137, #11018 and 
#12861) are ill-fit by either family of model SEDs ( X 2 > 5- 
11). This could result from systematic errors in our photometry 
for these galaxies (although we have checked that their colors 
do not appear to be strongly effected by nearby galaxies), or 
because their SEDs are not well-matched to the two simple al- 
ternatives used in our analysis: e.g. composite systems with 
evolved stellar populations and AGN contributions, galaxies 
with complex dust geometries, or very high-z starbursts which 
are also not included in our model (e.g. Waddington et al. 1999). 
We will discuss the properties of these sources in the context of 
a fuller photometric analysis of the complete radio population 
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in this field in Owen et al. (2002). 

Turning to the radio-undetected EROs we find that 20 of 
the 28 sources (79%) are reasonably described by one or other 
of the two families of model SEDs. Of these, 1 1 (39%) are 
better fit with a passive, evolved SED, compared to 9 (32%) 
which are well-fitted by a dusty SED. The median Zf-band 
magnitudes of the two classes are K = 19.60 ±0.11 for the 
evolved subsample and K = 19.67 ± 0.10 for the dusty sources, 
if we include with the latter the 1 1 well-fit dusty, radio-detected 
EROs (K = 19.91 ±0.28) we obtain a median magnitude of 
K = 19.82 ±0.09 for all of the EROs whose colors are well- 
fit by a dusty SED. Thus there is a slight tendency for the EROs 
with dusty SEDs to be fainter than those with evolved SEDs. 
The 8 radio undetected EROs which are not well-fit with either 
family of model SEDs (ignoring #2418 and #2549) are surpris- 
ingly bright on average K = 19.22 ± 0.24. To test whether this 
results from an error in the relative calibration of our photome- 
try we increase the estimate of the systematic errors in the rel- 
ative calibration of the photometry from Am = 0.05 to 0.2 and 
confirm that this does not significantly effect the average mag- 
nitude of the badly-fitting SEDs, although it does reduce their 
number from 8 to 4 - with the EROs making up the difference 
being equally classed as dusty and evolved systems. We also 
confirm that photometry of these badly-fit EROs is not unduly 
affected by crowding. 
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Fig. 5.— (R-K)-(J -K) color-color diagram for the 47 high-SNR 
ERO's brighter than K = 20.5 used in the SED fitting. We identify 
radio-detected and undetected EROs using filled and open symbols 
and in addition code the points depending upon the better fitting SED: 
either dusty, evolved or badly-fit. The solid line represents the bound- 
ary between dusty starburst and evolved, passive EROs proposed by 
Pozzetti & Mannucci (2000) - the dusty galaxies should lie to the 
right of the line (redder (J-K) colors), with the evolved systems 
on the left. Note that the high-SNR, radio-selected EROs lie either 
close to, or to the right, of this division, as do the majority of the 
high-SNR, radio-undetected EROs with dusty SEDs, confirming that 
the proposed boundary is reasonably effective. The dotted line di- 
vides low-mass stars from galaxies (which are redder in (J-K)). The 
dashed line shows the (R- K) = 5.3 boundary used to define our ERO 
sample. 

The median colors of the various subsamples of radio- 



undetected EROs are given in Table 2. The dispersion of the 
colors of the individual EROs in each subsample around their 
median are also listed there, these show that the dispersion 
within the subsamples of radio undetected EROs classed as 
dusty or evolved is less than that seen for the whole sample. 
This confirms that as expected the photometric classifications 
isolate more homogeneous subsamples from the overall radio- 
undetected ERO population, with the exception of the badly-fit 
subsample whose larger dispersion suggests that it still contains 
a mix of SED types. 

3.5.2. The mix of dusty and evolved EROs 

The strong preference for fitting the colors of high-SNR 
radio-detected EROs with the dusty SED model, 11/12 well- 
fit examples, gives us confidence that this model can be re- 
liably employed to identify similarly dusty, active systems in 
the radio-undetected ERO sample and hence estimate the preva- 
lence of such galaxies in the full ERO population. 

Taking the estimate that at least 20% (9/45) of the radio- 
undetected EROs have dusty SEDs and including the 1 1 high- 
SNR, radio-detected EROs with similar photometric classifica- 
tions, we conclude that at least > 44% (20/45) of those EROs 
with well-sampled SEDs have colors compatible with dusty 
starbursts, or > 30% (20/66) of all EROs with (R-K)> 5.3 
and K < 20.5 (median K <~ 19.6). This is consistent with the 
estimate of the proportion of dusty starbursts at K < 19.2 (me- 
dian K ~ 18.3) of 33-67% from Cimatti et al. (2002). 

There are two different assumptions which can now be made 
to broaden this constraint: one addresses the failure rate for fit- 
ting dusty EROs with our model SED, the other with the num- 
ber of dusty systems in the low-SNR ERO sample. Firstly, we 
could assume that the proportion of dusty SEDs in the high- 
SNR radio-detected sample is the same as that measured for 
those EROs with well-fitted SEDs, 92% (11/12), i.e. that the 
majority of the radio-detected EROs have dusty SEDs. Based 
on this we can then estimate the failure rate for fitting the col- 
ors of a dusty ERO with our model template as 5/16. Adopting 
the same failure rate for the dusty galaxies in the high-SNR 
radio-undetected sample, which have similar quality photom- 
etry, means that half of the badly-fit EROs are actually dusty 
galaxies (we would get the same estimate if we had simply pro- 
portioned the badly-fit ERO sample on the basis of the relative 
numbers of the successful fits). Hence we would estimate that at 
least 46% ( 1 3/28) of the high-SNR radio-undetected EROs have 
dusty SEDs. Adding in the radio-detected EROs then leads to 
a lower limit of > 50% on the proportion of dusty starbursts in 
the full ERO sample. Secondly, we could assume, based on 
the overall similarity of their median colors, that the high-SNR 
EROs are representative of the full samples. If this assumption 
is adopted then we would again conclude that a minimum of 
> 44% of the EROs in our full sample have colors which are 
well-fit by a dust-reddened starburst. 

Hence, either of these assumptions would suggest that a min- 
imum of about a half of the (R-K)> 5.3 and K < 20.5 ERO 
population have colors consistent with those expected for dusty 
starbursts. As these two assumptions are independent, we can 
also combine them by applying the mis-classification correc- 
tion to the radio-detected and undetected samples (giving 16/17 
and 13/28) and then assuming that the relative proportions of 
EROs of the photometric types in the high-SNR samples apply 
independently to the full samples (yielding 20/21 and 21/45). 
In this way we estimate that perhaps 60 ± 15% of the full ERO 
sample are distant, dusty starburst/AGNs, with the remainder 
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Fig. 6. — The broad-band SEDs for the high-SNR, radio-detected and undetected EROs from the sample used in our photometric analysis, each 
panel is ordered in terms of increasing (J-K) color. The left-hand and middle panels show EROs which are undetected in our radio map, while the 
right-hand panel shows the SEDs for the 17 high-SNR, radio-detected EROs. These SEDs are plotted on an arbitrary flux scale and we indicate 
uncertain regions of the SEDs by dashed lines. We also show an example of a dusty starburst SED and a passive, evolved galaxy, both at z = 1 
at the bottom of the central panel, as well as the central wavelengths of the 6 passbands used in our analysis. Note the peculiar SEDs of ERO 
#2418 and #2549, we identify #2418 as a late-type M dwarf (as suggested by its colors in Fig. 6), while the comparison of our independent J- and 
AT-band observations of the field shows that #2549 is a transient source. 



being passive systems (although these may also include a small 
number of very high-redshift galaxies, e.g. Waddington et al. 
1999). This estimate of the proportion of evolved and dusty 
EROs in the K < 20.5 population is of course also sensitive to 
the strong clustering of the passive ERO population. 

Our photometric classification also provides some other in- 
formation about the EROs: in particular the best-fit redshift and 
reddening. However, the limited precision of our photometry 
and the strong degeneracies between these parameters means 
that they should be viewed with caution. For the high-SNR, 
radio-detected EROs, which are well-fit by the dusty SED, the 
median redshift is z = 1.0 ±0.3, with a quartile range of z = 0.8- 
1.5. This range is very similar to that seen in the spectroscopic 
survey of Cimatti et al. (2002) and gives us some confidence in 
the average properties we determine for this population. The 
median reddening for this sample is Ay ~ 3.5 ± 0.4. For the 
radio-undetected EROs, those galaxies which are fitted with 
a dusty SED give a similar median redshift of z = 1.0 ±0.2 
and a slightly lower reddening of Ay <~ 2.0 ±0.5. These es- 
timated reddenings are consistent with that estimated for their 
composite spectrum of star-forming EROs by Cimatti et al. 
(2002), Ay ~ 2.4 ± 0.9. For both subsamples the median 
TT-band luminosities at the estimated redshifts corresponds to 
My ~ -18.4 ± 0.4, or My <~ -21 .5 when corrected for redden- 
ing. The unobscured luminosities of these galaxies would be 
comparable to L v galaxies today, however, their young stellar 
populations would be expected to fade substantially and so it is 
clear that these galaxies are likely to evolve into sub-L v galax- 
ies at the present-day, unless the current starburst continues for 
> 10 8 yrs. 

Next, we look at the relationship between our crude restframe 
visual morphologies and the preferred photometric classifica- 
tion for 28 EROs (8 radio-detected and 20 radio-undetected, 
Table 1). The morphological mix in both cases is similar 



to the proportions of the different SED types, if we inter- 
pret the compact EROs as passive, evolved systems (Man- 
nucci et al. 2002). For the radio-undetected sample we find 
a broad equality between the classes: 29%:29%:42% for com- 
pact:disky/amorphous:faint, similar to 39%:32%:29% for the 
evolved, dusty and badly-fit SEDs. Similarly, for the radio- 
detected sample, which is dominated by dusty SEDs, the mor- 
phological mix is: 0%:82%:18% compared to the distribution 
of SED fits: 6%:65%:29%. However, a one-to-one comparison 
for individual EROs shows that while the distributions appear to 
be telling the same story, the individual galaxies are much more 
mixed: half of those EROs with dusty SEDs have disky, amor- 
phous or disturbed morphologies, but none of the EROs with 
passive SEDs have compact morphologies (this morphological 
class is almost evenly split between dusty SEDs and failed fits). 
We also note that there is no clear trend in the SEDs of those 
radio-detected EROs which are resolved in our VLA map. 

To summarize, we place a firm lower limit of > 30% on 
the proportion of dusty starbursts in the full (R-K) > 5.3 and 
K < 20.5 ERO population. Adopting two separate corrections 
for the number of dusty EROs which are missed either due to 
the failure of our photometric classification or because they 
were removed due to their poor photometry, we propose that 
> 45-50% is a more realistic lower limit on the fraction of dusty 
EROs. Finally, by combining these corrections we suggest that 
the actual fraction of EROs with dusty starburst/ AGN-like col- 
ors is probably close to 60 ± 15%. 

4. DISCUSSION 

4. 1 . The evolution of dusty and evolved EROs 

We have estimated the relative proportion of dusty, active and 
passive, evolved systems in the K < 20.5 ERO population. By 
comparing this to the mix of classes estimated for the K < 19.2 
sample by Cimatti et al. (2002) we can identify any decline in 
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either class of EROs which might be responsible for the ap- 
parent break in the ERO counts at K ~ 19-20 (McCarthy et 
al. 2001; Smith et al. 2002a). The similarity of the two esti- 
mates suggests that there is no evidence for a significant shift 
in the mix of these two sub-populations across K ~ 18-20, al- 
though it should be noted that the Cimatti et al. sample is se- 
lected using a slightly bluer color criteria, (R-K) > 5, which 
complicates this comparison. This apparent lack of a substan- 
tial shift in the nature of the ERO population may indicate that 
the break arises instead from the volume sampled - with K < 20 
EROs populating a relative narrow redshift range. This volume 
would be bounded at the low redshift end by the requirement 
that the galaxy is observed at high enough redshift, i.e. suf- 
ficiently far into the restframe UV, to enable a dust-reddened 
starburst or evolved stellar population to appear extremely red, 
(R-K) > 5.3 (Daddi et al. 2001; Firth et al. 2002). Observa- 
tionally this appears to require z Z 0.8-1 (§3.5; Cimatti et al. 
2002). At the upper-end the redshift limit may be imposed ei- 
ther by the luminosity function, passive bluing of the stellar 
population, or the increasing fraction of unobscured star forma- 
tion in both sub-classes of ERO (the continuity between high-z 
SCUBA galaxies and the radio-selected ERO population sug- 
gested below would argue against the latter suggestion). If the 
break in the counts is due to the surveys reaching beyond the 
volume containing the ERO population then the K ~ 19 break 
(Fig. 2) should roughly correspond to M* for a composite lumi- 
nosity function with a Schechter form. This would imply that 
the characteristic luminosity of the ERO population is very sim- 
ilar to that seen in the local field, M v ~ -21.5, and that the ERO 
luminosity function has a faint end slope, a, which is shallower 
than a ~ -0.7 (the slope of the differential counts beyond the 
break). To illustrate this we plot a toy model in Fig. 2, this 
has a mix of Gaussian luminosity functions each with mean 
luminosity of M K = 24.3 (M* +0.7 at z = 0, Cole et al. 2001) 
and a dispersion of a K = 1 (similar to the width of the lumi- 
nosity function of elliptical galaxies locally). These are com- 
bined to produce a redshift distribution with a mean redshift of 
<z>^ 1.1 and a dispersion of a z = 0.2 (Cimatti et al. 2002; 
Firth et al. 2001), and the cumulative counts calculated using a 
simple K-correction. Clearly this toy model (with some further 
tuning) can reasonably reproduce the ERO counts, indicating 
that a combination of a narrow luminosity function and a nar- 
row redshift range are a sufficient, if not unique, explanation 
of the form of the counts of the whole ERO population across 
K~ 17-21. 

However, such a model has two obvious shortcomings - 
firstly, it requires that the evolutionary behavior of the dusty 
and passive ERO populations is similar, which seems highly un- 
likely, and secondly the required evolution means that the dusty, 
active EROs selected in our radio map are distinct from the 
higher redshift and fainter counterparts to the SCUBA galax- 
ies, rather than representing a simple continuity of the popula- 
tion to lower redshifts and brighter K-banA fluxes. For these 
reasons we reject this model and construct a second toy model. 

For the second model we assume that the number counts 
of dusty EROs are well-described by the power-law slope es- 
timated for the radio detected population in §3.1, but with a 
normalization increased to reproduce the ~ 60% fraction of 
the K < 20.5 ERO population which we estimate have dusty 
SEDs. To this we add a Gaussian population representing the 
passive sub-population and normalize the combined sample to 
reproduce the weighted cumulative surface density of EROs at 
K < 20.5 from Thompson et al. (1999) and Smith et al. (2002a). 



Representative parameters for the passive sub-population re- 
quire a mean magnitude of K = 19.1 (M£ + 0.5 at z = assum- 
ing no evolution, or M* K +1.5 assuming passive evolution from 
z ~ 1.1) and a FWHM of 1.2 mag. We illustrate the cumulative 
counts of both sub-populations and the combined population in 
Fig. 2. We note that this model predicts that the relative pro- 
portions of dusty and evolved EROs for a K < 19.2 sample is 
35%:65%, in agreement with the current limits from Cimatti et 
al. (2002). Clearly this model, even without elaborate tuning, 
can provide an adequate description of the current data, while 
at the same time providing continuity between the dusty EROs 
detected in our VLA map and the higher redshift and fainter 
counterparts to submm sources (Smail et al. 2002; Ivison et al. 
2002) and we therefore prefer it. A clear distinction between 
the alternative models will come from extending the counts of 
EROs fainter, to K > 21-23, and estimating the relative mix of 
evolved and dusty systems at these depths. However, the depth 
of the imaging required for such a survey, R <~ 28 and K ~ 23, 
will make it observationally demanding. 

4.2. Star-formation in the dusty ERO population 

The resolution of over a third of the radio-detected EROs 
in our VLA map, along with the close correlation between the 
radio and optical morphologies of these EROs, suggests a com- 
mon mechanism for the emission at these two wavelengths. The 
simplest explanation is that the emission in both bands is pow- 
ered by massive star formation within an extended region - on a 
scale of ~ lOkpc. Thus for at least a third of the radio-detected 
EROs, star formation is likely to contribute the majority of the 
radio emission (this includes one of the brightest radio EROs, 
#01053). Furthermore, we note that using the median 1.4-GHz 
surface brightnesses of these resolved EROs and the star forma- 
tion calibration discussed below, they are expected to have high 
star formation surface densities: ~ 4M Q yr" 1 kpc" 2 , more than 
sufficient to drive powerful winds. 

4.2. 1 . The fraction of AGN- dominated EROs 

To reliably quantify the star formation occurring in the com- 
plete radio-detected ERO sample we must estimate the contam- 
ination from AGN-dominated systems. The AGN fraction in 
the whole /iJy radio population can be determined from the 
spectroscopic and morphological surveys of fiJy radio sources 
(e.g. Muxlow et al. 1999; Roche, Lowenthal & Koo 2002). 
These suggest that at most 20% of radio sources at > lOO^Jy are 
likely to be AGN. However, these estimates are typically based 
on samples of radio sources with brighter optical counterparts 
than the ERO sample analysed here (and the AGN tend to be 
the bluer sources in even these samples, Roche, Lowenthal & 
Koo 2002; Waddington et al. 2000). Nevertheless, non-thermal 
emission is seen in a handful of well-studied EROs (Pierre et al. 
2001; Willott et al. 2001; Smith et al. 2001; Cowie et al. 2001) 
indicating that AGN are not completely absent from the ERO 
population, but the fraction of these in complete photometri- 
cally selected ERO samples is not well determined (Cimatti et 
al. 2002). The most useful constraint on this fraction for our 
purpose is that estimated from the number of ERO counterparts 
to faint X-ray sources (Alexander et al. 2001, 2002; Brusa et 
al. 2002; Mainieri et al. 2002). While we earlier confirmed that 
none of the EROs in our sample is associated with a luminous 
X-ray source, using archival XMM-Newton observations, much 
more stringent limits are now available on X-ray counterparts 
to EROs from recent Chandra and XMM-Newton observations 
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of blank fields (which have lower backgrounds than our cluster 
field). 

Alexander et al. (2002) identify six (I-K) > 4 and K < 20.1 
X-ray sources in a 70.3 sq. arcmin region within the 1-Ms 
Chandra Deep Field North exposure, giving a surface density 
of 0.085 ± 0.035 arcmin" 2 . Based on the relative numbers of 
(I-K) > 4 and (R-K) > 5.3 sources and the number counts of 
(R-K) > 5.3 EROs we would expect roughly - 0.07 ±0.03 
arcmin" 2 X-ray-detected EROs with (R-K)> 5.3 and K < 
20.5. This estimate assumes that the X-ray-detected fraction 
is not a strong function of color or ^-magnitude, although we 
note that 2/6 sources in the Alexander et al. (2002) sample have 
(I-K) = 4.0 and hence are borderline cases. Compared to our 
ERO surface density of (1.04±0.12) arcmin" 2 this indicates 
that around ~ 10% of the ERO population are detected in very 
deep X-ray observations. 

However, the Chandra exposure used by Alexander et al. 
(2002) is sufficiently sensitive that it will detect not only ob- 
scured AGN sources, Lx Z 10 44 ergs s" 1 , but also powerful star 
forming galaxies at z ~ 1 (Ranalli et al. 2002). If we in- 
stead restrict the X-ray sample to those EROs which are too 
X-ray luminous to be ULIRG-like starbursts at z ~ 1, i.e. 
Lx(0.1-2.4KeV) > 4 x 10 42 ergs s" 1 (twice the X-ray luminos- 
ity of the local X-ray-bright starburst NGC 3256, which has 
log 10 (L F/R /L©) ~ 11.8, Moran et al. 1999), then the number 
of probable AGN-dominated sources in the Alexander et al. 
(2002) sample would drop to three (two of which have bor- 
derline, (I-K) = 4.0, colors), or an estimated surface density 
of < 0.03 arcmin" 2 using our ERO criteria. We conclude that 
at most 3% of the ERO population with (R-K) > 5.3 and 
K < 20.5 are likely to host energetically dominant AGN. This 
is consistent with the estimates from ERO counterparts to hard 
X-ray sources detected by XMM-Newton in Franceschini et al. 
(2002), as well as the recent estimate that 1-2% of the ERO 
population have hard X-ray counterparts by Brusa et al. (2002) 
and Mainieri et al. (2002). 

While the occurrence of AGN-dominated sources in the ERO 
population is obviously low, the estimates derived above do not 
tell us how frequently radio-selected EROs are detected in the 
X-ray waveband. Fortunately, we can also crudely estimate this 
from the Alexander et al. (2002) survey. Three of the X-ray- 
detected sources with (I-K) > 4 in Alexander et al. (2002) are 
also detected in deep 1 .4-GHz observations, two of these have 
X-ray and radio emission consistent with an X-ray bright star- 
burst at z ~ 1, suggesting a surface density of ~ 0.02 arcmin" 2 
for X-ray-detected AGNs with /iJy ERO hosts, corresponding 
to ~ 6% of the whole radio-detected ERO population. This rate 
of AGN contamination is consistent with our identification of 
only one radio-detected ERO with an evolved SED (assuming 
the radio emission from this galaxy arises from a weak AGN, 
see §3.5), corresponding to 6% of the sample. We therefore 
adopt 6% as the AGN contamination in our sample and conser- 
vatively assume that the radio emission from any EROs hosting 
an AGN is dominated by the AGN. 

4.2.2. The star formation density in EROs 

We can now use our ERO sample to estimate the star forma- 
tion density in the most obscured starbursts at z ~ 1 (Cram et 
al. 1998; Oliver et al. 1998; Haarsma et al. 2000). We calculate 
two estimates of the star formation density - first we take the 
minimum number of galaxies in our sample which are likely to 
be dusty, active systems based on the analysis of their optical- 
near-infrared SEDs: this subsample comprises 1 1 high-SNR, 



radio-detected EROs and 9 radio non-detections, and provides 
a conservative lower bound on the estimated star formation den- 
sity. 

Those radio-undetected EROs with dusty SEDs are likely to 
either be more distant or have typically lower star formation 
rates (and hence radio fluxes) than the radio-detected examples. 
Based on the analysis in §3.5.2 there is a hint that these galax- 
ies have similar redshifts to the radio-detected EROs, but lower 
reddening's - and hence perhaps lower star formation rates 
(Hopkins et al. 2001). Using their estimated surface density and 
extrapolating the cumulative number counts of radio-detected 
EROs in our field from §3.2, we would expect most of the dusty, 
radio-undetected EROs to have radio fluxes above ~ 5^Jy and 
hence we assume that the typical dusty, radio-undetected ERO 
has a radio flux of ~ 10 /zJy. 

For the second estimate we adopt the incompleteness esti- 
mates for the photometric classification from §3.5 and hence 
use 20 radio detections and 21 radio-undetected EROs. Sim- 
ply summing the 1 .4-GHz fluxes for these two subsamples we 
obtain a total radio flux of > 910 fiJy for the conservative case 
and 1370 [iJy for the more realistic estimate (these include con- 
tributions from the radio-undetected starbursts of 90^iJy and 
210 /iJy respectively). Applying our estimated AGN contami- 
nation reduces both of these estimates by 6%, however, if this 
contamination preferentially effects the most luminous, unre- 
solved radio EROs then the radio fluxes could be reduced by as 
much as <~ 15% (roughly half the difference between our "con- 
servative" and "realistic" cases). 

Look— back time (Gyr) 
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Fig. 7. — The evolution of the star formation density in the Universe 
as measured using a range of long-wavelength tracers (which should 
be relatively insensitive to dust obscuration). We plot our measure- 
ment of the star formation density in the dusty ERO population, the 
lower limit of the error-bar on this point corresponds to our conser- 
vative estimate of the star formation density in this population (see 
text). We compare this to other estimates from the radio (Haarsma et 
al. 2000; Serjeant et al. 2000); far-infrared/submillimeter (Blain et al. 
1999a); mid-infrared (Flores et al. 1999); near-infrared (Connolly et 
al. 1997), Ha (Gallego et al. 1997; Yan et al. 1999) and UV/optical 
(Lilly et al. 1995). 

Hence we estimate cumulative restframe 1 .4-GHz luminosi- 
ties from star formation in the two subsamples of L\ 4 > (8.0 ± 
1.3) x 10 24 WHz and L XA = (12.2 ± 1.5) x 10 24 WHz respec- 
tively. These assume a K-correction based on a v~ Q1 spectrum 
typical of star forming galaxies (Richards 1998), which gives 
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a correction of 1.62^2, where the error is dominated by the 
redshift range of the galaxies. We now use the local correla- 
tion between 1.4-GHz luminosity and far-infrared luminosity 
(Condon et al. 1991) to convert this estimate into a cumulative 
far-infrared luminosity for the sample, we note that the validity 
of this correlation has recently been confirmed at the redshifts 
studied here by Garrett (2002). In this way we obtain cumula- 
tive far-infrared luminosities of: L FIR > (16.1 ±2.5) x 1O 12 L 
and (24.3 ±3.7) x 1O 12 L respectively. We estimate that these 
correspond to SFRs forM > 5M of > 1450 ± 220 M yr" 1 and 
2180 ± 33OM yr" 1 using the conversion in Condon (1992). If 
we instead adopt the new calibration of Yun, Reddy & Condon 
(2001), we would estimate star formation rates roughly half of 
the quoted values (see also Carilli et al. 2002). However, for 
comparison with previous work we retain the Condon (1992) 
calibration. Assuming a Salpeter IMF from 0. 1-100 M these 
rates are ~ 5.5 x higher. 

To calculate the volume probed by our ERO sample we use 
the photometric redshift estimates for the 20 radio-detected and 
undetected EROs which are classed as dusty starbursts and find 
a quartile range of z = 0.8-1 .5 (in reasonable agreement with the 
spectroscopically determined redshift range from Cimatti et al. 
2002, z = 0.8-1.4, see Daddi et al. 2002), which gives a corre- 
sponding volume of(l.l±0.2)x 10 5 Mpc 3 at a median redshift 
of z = 1 • The quoted uncertainty corresponds to changing either 
of the bounds on the volume by 5z = ±0.2 and encompases the 
estimate we would derive if we had instead adopted the redshift 
range for EROs from Daddi et al. (2002). 

Combining this volume with the star formation rates we esti- 
mate equivalent star formation densities of p# (O.1-1OOM ) > 
0.073 ±0.018M Q yr" 1 Mpc" 3 and 0.1 1 ± 0.03 M yr" 1 Mpc" 3 
(if the AGN contamination is predominantly in the most lumi- 
nous, unresolved radio sources the latter estimate will fall to 
~ 0.09 M yr" 1 Mpc" 3 ). We compare these measurements to 
other estimates of the star formation density in Fig. 7. 

Our measurement of the star formation density in the radio- 
detected ERO population at z ~ 1 is comparable to estimates of 
that in Ha-selected galaxies at these epochs (Yan et al. 1999), 
and greater than that seen in optically-selected galaxies (Lilly 
et al. 1995; Connolly et al. 1997). There is little or no over- 
lap between the ERO population used in our analysis and the 
classes of galaxies which are employed in these previous stud- 
ies and hence their contributions to the total star formation den- 
sity are distinct. This implies that the activity in obscured galax- 
ies makes a sizable contribution to the total star formation at 
Z ~ 1, especially when allowance is made for the proportion of 
the obscured star forming population which are missing from 
our ERO sample due to the escape of a small fraction of UV 
light. The significant contribution from obscured star forma- 
tion at z ~ 1 confirms earlier claims by Flores et al. (1999) and 
Haarsma et al. (2001) and is particularly important as this era 
represents the peak in the star formation activity in the Universe 
(note that when integrated over cosmic time, the distribution in 
Fig. 7 is highly peaked at z ~ 1-2). This result therefore ex- 
tends the conclusions of the mid-infrared, submillimeter and ra- 
dio surveys, that the use of obscuration-independent measures 
of the star formation rate is critical to reliably estimate the total 
star formation density, from z ~ 2-3 to the era when most of 
the stars seen in galaxies today were formed. 

4.3. Dusty EROs in the Far-infrared and Submillimeter 

We now estimate the volume density of the most strongly 
star forming EROs. Based on the photometric redshifts 



for the 11 radio-detected EROs which are well-fit by a 
dusty SED, we estimate a median radio luminosity of Z4.4 = 
10 23 ' 6±0 ' 3 WHz" 1 (compared to the Galaxy's radio luminosity 
of Z4.4 = 10 21 4 W Hz" 1 ). Assuming these galaxies follow the lo- 
cal radio-far-infrared correlation, this translates into a median 
far-infrared luminosity of Lfir/Lq = I0 n ' 9±a3 for the subsam- 
ple. The brightest five EROs in this subsample have inferred 
far-infrared luminosities above L FIR > 1O 12 L which classes 
them as ultraluminous infrared galaxies (ULIRGs). Based on 
these five galaxies we estimate a surface density of radio- 
detected ULIRG EROs of 0.07 ±0.03 arcmin" 2 . This estimate 
could be increased by 1 .4 x to account for radio -detected EROs 
which are not well-fit by our model SEDs, although equally any 
increases should be offset against any AGN contribution to the 
radio luminosities of these galaxies. 

Using the volume containing the dusty ERO sample given in 
§4.2 we estimate an equivalent space density of ULIRG EROs 
as - (0.5 ±0.2) x lO^Mpc" 3 (and - (2.7 ±0.4) x lO^Mpc" 3 
for all dusty EROs brighter than K < 20.5). This is significantly 
higher than the space density of comparable luminosity events 
locally, 0.03 x 10" 4 Mpc" 3 (Sanders & Mirabel 1996), even be- 
fore allowance is made for the fact that some high-z ULIRGs 
will be missing from our ERO selection due to their blue colors 
(Trentham et al. 1999; Goldader et al. 2002). Based upon the 
far-UV colors of local ULIRGs it appears that more than half 
of this population would be too blue at z ~ 1-2 to be included 
in our ERO sample. Hence, if local ULIRGs are good analogs 
for the events which produce the high-z ULIRG population then 
our volume density of z ~ 1 ULIRGs would require an upwards 
revision by a factor of > 2x. We conclude that there has been 
at least an order of magnitude decrease in the number density 
of these luminous infrared events between z ~ 1-2 (8-10 Gyrs 
ago) and the present-day. This trend is consistent with the evo- 
lutionary models of Blain et al. (1999a), which also fit the num- 
ber counts of SCUBA galaxies at higher redshifts (Fig. 7). 

We can also estimate the submillimeter fluxes from the radio- 
detected EROs and hence investigate their relationship to the 
850-/im selected samples detected by SCUBA. Based on the 
individual radio fluxes and estimated redshifts for the EROs we 
would predict a median 850-/im flux of 2 mJy and a range of 
0.3-5 mJy (Carilli & Yun 2000). This places the bulk of these 
galaxies below the flux limit of typical SCUBA observations of 
faint radio sources (e.g. Chapman et al. 2001) and at least half 
of them below the confusion limit of the deepest (unlensed) 
SCUBA surveys. Clearly the radio/ERO selection allows us to 
probe the star forming population to lower star formation rates 
(and also lower redshifts) than typically achieved in submm sur- 
veys as demonstrated by our detection of all three cataloged 
submm sources in our field (see Ledlow et al. 2002). 

Nevertheless, the most luminous radio-detected EROs should 
overlap with samples of ERO counterparts to SCUBA galaxies. 
We can test this by comparing the surface densities of these 
two populations. Our estimate of the surface density of > 2 mJy 
radio-detected EROs is 0.07 ± 0.03 arcmin" 2 at K < 20.5. This 
can be compared to the claim that 13 ±9% of the counterparts 
to SCUBA galaxies brighter than 2mJy have K < 20.5 and 
(R-K) > 5.3, equivalent to a surface density of ~ 0.11 ±0.09 
arcmin" 2 , (Smail et al. 1999a, 2002; see also Mohan et al. 2002) 
and the estimate surface density of EROs brighter than 8 mJy of 
0.03 arcmin" 2 from Ivison et al. (2002). Allowing for the large 
uncertainties in these estimates, and the possible upward re- 
vision of the density of radio-selected ULIRGs to account for 
badly-fit SEDs (including the known SCUBA ERO #01206), 
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we conclude that they are in reasonable agreement. This sug- 
gests that by combining sensitive radio and submillimeter sur- 
veys it will be possible to track the evolution of the star forma- 
tion density in highly-obscured galaxies from z = to z > 3-5 
(Owen et al. 2002; Chapman et al. 2002b). 

Finally, we estimate the contribution from the dusty ERO 
population to the far-infrared background at a wavelength cor- 
responding to the peak of their dust emission: A = (1 +z)X pe ak ~ 
120/im. Scaling from the 1.4-GHz emission we estimate that 
the dusty ERO population produces vl v ~ 1 x 10~ 9 Wm~ 2 sr" 1 
at <~ 100/im, compared to estimates of the background at this 
wavelength of vl v = (2 ± 1) x 10" 8 W m 2 sr 1 (Hauser et al. 
1998; Lagache & Puget 2000). Thus the dusty EROs could 
contribute <~ 5-10% of the far-infrared background at its peak at 
~ 100/im. This compares to the SCUBA population which pro- 
duces - 100% of the background at 850/mi (Blain et al. 1999b; 
Smail et al. 2002; Cowie et al. 2002) and probably the bulk of 
the far-infrared background at A > 200/mi (where the dust emis- 
sion peaks for a population with a characteristic temperature of 
40 K at z ~ 2-3). Both the dusty ERO and the SCUBA popula- 
tion contribute <C 1% of the extragalactic background at 1/xm, 
underlining the very different nature of the populations selected 
in the optical and the submillimeter/radio wavebands. 

4.4. The Passive ERO Population 

Our suggestion that at least a half of the faint ERO pop- 
ulation are dust-enshrouded starbursts obviously reduces the 
estimated volume density of passive, evolved EROs. We es- 
timate that passive EROs at z > 1 have a volume density of 
~ 3.6 ±0.6 x 10" 4 Mpc" 3 for galaxies with a typical luminosity 
of My ~ -20.6 ± 0.4. Passive evolution of the stellar popula- 
tions would produce ~ 0.9 mags of brightening in the observed 
TT-band to this redshift, meaning that these galaxies would have 
M v ~ -19.7 ±0.4 at the present day if they evolve passively. 
This assumes that the passive EROs inhabit the same volume 
estimated for the dusty ERO population, although as we dis- 
cussed in §4. 1 these systems may lie in a smaller volume and 
hence this volume density should be taken as a lower limit. 
In this regard we note that the K < 19.2 passive sample from 
Cimatti et al. (2002) have z ~ 1.04 ±0.15, where the error 
represents the l-cr dispersion, and a quartile range of z = 0.8- 
1.2, somewhat narrower than the redshift range we use - and 
suggesting a volume for the sample closer to 0.6 x 10 5 Mpc 3 . 
Adopting this volume would roughly double our estimate of 
the volume density of passive EROs. 

The most frequently proposed local analogs to passive EROs 
are elliptical galaxies, although it is clear that apparently pas- 
sive EROs embrace a wider cross-section of the z ~ 1 popula- 
tion than just disk-less bulges (Smith et al. 2002b). It is con- 
sequently difficult to compare the volume density of this pop- 
ulation with that of their possible descendants at the present- 
day. However, we note that the combination of our reduction 
in the proportion of passive galaxies in the ERO population, 
coupled with the increase in the diversity of their morphologi- 
cal counterparts identified by Smith et al. (2002b), means that 
previously well-fitting pure luminosity evolution models with 
a high redshift of formation, z > 2, will now over predict the 
ERO counts at K ~ 20 (Daddi et al. 2000; Firth et al. 2001; 
Smith et al. 2002a). This suggests there may be some scope 
for the recent formation (or transformation) of some fraction of 
the local elliptical galaxy population (e.g. Aguerri & Trujillo 
2002). Nevertheless, the reduction in the proportion of passive 
systems in the ERO population is still not sufficient to recon- 



cile the observed surface density of evolved EROs with that 
predicted by semi-analytic galaxy formation models (Smith et 
al. 2002a; Cole et al. 2000). 

The proposed reduction in the passive ERO fraction also has 
ramifications for several previous results based on the proper- 
ties of the whole ERO population. The most obvious is the in- 
terpretation of the strong clustering signal seen in ERO samples 
(Daddi et al. 2001, 2002; McCarthy et al. 2001). The strong 
angular clustering of the ERO population, especially for the 
brightest examples, appears to vindicate the predictions of hier- 
archical galaxy formation models assuming these galaxies rep- 
resent the passive, evolved progenitors of present-day ellipticals 
(Firth et al. 2001). In the same manner, the apparent agreement 
between the angular clustering of Lyman-break objects at z ~ 3 
and theoretical models of halo clustering has provided support 
for the suggestion that these objects are the progenitors of mas- 
sive ellipticals (Giavalisco & Dickinson 2001). 

However, the presence of a large proportion of, perhaps 
weakly clustered, dusty EROs in the bright ERO population 
would tend to dilute the angular clustering signal arising from 
the strongly clustered passive EROs, unless there is a strong 
cross-correlation between the two populations (see Daddi et al. 
2002 and §3.3). Indeed, a further element comes into play when 
attempting to interpret the projected clustering of ERO sam- 
ples: the redshift depth of the volume containing the galaxies. 
If the passive ERO population lies in a relatively narrow redshift 
range (see the toy models in §4.1; Firth et al. 2001) then the di- 
lution effect from the depth of the sample will be less than an- 
ticipated, requiring an intrinsically less strongly clustered pop- 
ulation to produced a given angular clustering signal (Daddi 
et al. 2001). Thus the effects of a narrower redshift depth for 
the evolved ERO subsample and the increased proportion of 
weakly-clustered, but typically faint, dusty EROs may cancel 
each other out leaving the estimated correlation length for the 
bright, evolved EROs close to that required for the progenitors 
of the L* elliptical galaxies seen in the local Universe, as esti- 
mated by Daddi et al. (2001) and McCarthy et al. (2001). We 
conclude that the complex interplay between the sample vol- 
ume and the mix of evolved and obscured galaxies in the ERO 
population as a function of apparent magnitude means that the 
interpretation of the angular clustering of this class of galaxies 
is far from trivial. 

5. CONCLUSIONS 

The main conclusions of this work are: 

1) We catalog 68 EROs within a 72.8 arcmin 2 field which has 
the deepest 1 .4-GHz observations of any region on the sky. We 
find that 31% of these EROs are detected in the radio map 
down to a flux limit of 12.6 /iJy, of these radio-detected EROs, 
33% are resolved in our VLA map with a typical FWHM of 
1 .4 ± 0.4" - equivalent to 12 ± 3 kpc at z ~ 1 . Moreover, the ra- 
dio emission in these EROs mirrors their /?-band morphologies, 
suggesting that star formation is responsible for the bulk of the 
radio flux. 

2) Using our deep RIzJHK imaging we photometrically clas- 
sify the radio-detected and undetected EROs into dusty or 
evolved systems based on simple SED models. We first analyse 
the colors of those radio-detected with high quality photometry 
and conclude that the majority of these are adequately described 
by a highly dust-reddened starburst spectrum. This provides 
support for applying this analysis to the radio-undetected EROs 
to search for similar obscured, active galaxies. 

3) Using our simple photometric models we estimate that more 
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than 32% of the radio-undetected EROs with high quality pho- 
tometry, or a minimum of 20% of all radio-undetected EROs, 
have the colors of dusty starbursts, similar to the radio-detected 
EROs. We estimate the rate of mis-classification of EROs in our 
analysis assuming that most of the radio-detected EROs are ac- 
tually dusty starbursts, and based on this propose that <~ 45% of 
the EROs with radio fluxes < 12.6/iJy have dusty SEDs. From 
the continuity of the number counts of radio-detected EROs we 
would expect that most of these sources would have 1 .4-GHz 
fluxes around 5 p]y. We also find that one of the 68 EROs has 
an SED consistent with galactic star, while another apparent 
ERO appears to be a transient source. 

4) Combining the samples of radio-detected and undetected 
EROs which are photometrically classified as dusty we place a 
firm lower limit of > 30% on the proportion of dusty starbursts 
in the whole ERO population. To obtain a more realistic esti- 
mate we adopt a correction for mis-classification and assume 
the resulting proportions of dusty and evolved EROs are repre- 
sentative of the full sample (not just those with the best pho- 
tometry) we estimate that at least 45-50%, and possibly up to 
~ 60 ± 15%, of the whole ERO population with (R-K) > 5.3 
and K < 20.5 are probably dusty, active galaxies. Our fitting 
procedure also suggests a median redshift of z ~ 1 .0 ± 0.3 with 
a quartile range of z = 0.8-1.5. This is in reasonable agree- 
ment with the spectroscopic survey of somewhat brighter and 
bluer EROs by Cimatti et al. (2002). Assuming that AGN do 
not dominate the radio emission from these galaxies, we esti- 
mate a volume density of EROs with far-infrared luminosities 
of > 1O 12 L (i.e. ULIRGs) of 0.5 x lO^Mpc" 3 , this is over an 
order of magnitude higher than for equivalent luminosity events 
at z = and confirms the strong evolution of the most dusty, ac- 
tive galaxies out to z > 1 (Blain et al. 1999a). 

5) We use the local far-infrared-radio correlation to estimate 
star formation rates from the radio fluxes of the dusty EROs, as- 
suming that their radio emission is purely powered by massive 
star formation. Combining these estimates with the probable 
volume spanned by our dusty ERO subsample we can calculate 
the star formation density in this obscured population. Adopt- 
ing conservative assumptions we find that the star formation 
density in star forming EROs at z ~ 1 probably exceeds p\ (0. 1- 
100 M Q ) > 0.073 ±0.018 M yr" 1 Mpc" 3 and maybe as high as 
p\= 0.1 1 ± 0.03 M Q yr" 1 Mpc" 3 . This measurement is compa- 
rable to estimates of the star formation density in Ha-selected 
galaxies at these epochs, and greater than that seen in optically- 
selected galaxies, showing that the activity in obscured systems 
may make a significant contribution to the total star formation 
at z ~ 1 . We also estimate that the dusty ERO population may 
produce up to 5-10% of the far-infrared background at its peak 
at ~ lOOjum, underlining the cosmological significance of these 
events for understanding the history of star formation in the 
Universe. 

6) We propose that the apparent break in the integrated counts 
of the EROs at K ~ 19-20 reflects the peak of the contribution 



from passive, evolved galaxies to this population. At magni- 
tudes fainter than K <~ 20 we predict that the population will be 
increasing dominated by dusty, active galaxies. The continuity 
in the properties of these faint dusty EROs connects them to 
samples of submm-selected obscured starbursts at higher red- 
shifts and fainter Zf-band magnitudes. This proposal can be 
tested by studies of the multiwavelength properties of extremely 
faint samples of EROs, K ~ 23, and their relationship to the 
blank-field populations of radio sources. 
7) By identifying the likely contribution from dusty EROs to 
the total ERO population we estimate that the volume den- 
sity of passive, evolved galaxies brighter than K < 20.5 is 
> 4 x 10" 4 Mpc" 3 . Recent results point to a morphological di- 
versity in the passive ERO population which would suggest that 
these galaxies should not be simply viewed as passively evolv- 
ing ellipticals (Smith et al. 2002b). We conclude that previ- 
ously well-fitting pure luminosity evolution models, which de- 
scribed the number counts of the whole ERO sample as due to 
a passively evolving elliptical population formed at high red- 
shifts, will probably over predict the number of truly passive 
EROs and hence there is scope for the recent (trans)formation 
of some local elliptical galaxies. Nevertheless, the presence 
of luminous, passive galaxies at z ~ 1 does point to a phase 
of massive spheroid formation at higher redshifts, consistent 
with the interpretation of the SCUBA population as the forma- 
tion phase of massive, proto-ellipticals (Lilly et al. 1999; Smail 
et al. 2002, 2003). The typical, reddening-corrected luminosi- 
ties of the radio-detected ERO population at z ~ 1 suggests that 
these galaxies will evolve into sub-L* systems by the present- 
day unless their star formation activity continues at a high rate 
for> 10 8 yrs. 

The area analysed in this work covers a small fraction of 
the sensitive field of our VLA map - we expect to be able to 
quadruple the ERO sample available for analysis in the near fu- 
ture. In addition we are undertaking a full photometric analysis 
of the <~ 1500 optically-bright p}y radio sources in this field 
using our panoramic UBVRIzJHK imaging dataset (Owen et al. 
2002), constrained by deep spectroscopy from GMOS on Gem- 
ini. 
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47 02 05.63 


20.11 


5.40 


3942 


09 42 27.216 


47 02 09.18 


19.31 


5.79 



37.5 
44.4 
39.4 
195.6 
112.2 
70.2 
51.9 
45.7 
44.1 
16.0 
24.4 
34.4 
64.5 
32.9 
81.1 
199.9 
23.2 
17.0 
24.6 
14.1 
16.3 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 

< 12.6 



a) Photometric errors in 0.01 mags are listed in []. b) Best-fitting model SED: D, dusty; E, evolved; B, badly fit. c) Morphological classification in K- and /?-bands: "C", compact/regular; 
"I", disky, disturbed or amorphous; "F", too faint to classify, d) SMM 109429+4658 (Smail et al. 1999). e) ERO #333 (Smail et al. 1999). /) ERO CI 0939+47 13 A (Persson et al. 1993). g) 
M star, h) Transient source. 



